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Abstract 

We propose a first-principles scheme for the description of the magneto-optical kerr effect within 
density functional theory (DFT). Though the computation of Kerr parameters is often done within 
DFT, starting from the conductivity or the dielectric tensor, there is no formal justification to this 
choice. 

As a first steps, using as reference materials iron, cobalt and nickel we show that pseudo- 
potential based calculations give accurate predictions. Then we derive a formal expression for the 
full dielectric tensor in terms of the density-density correlation function. The derived equation 
is exact in systems with an electronic gap, with the possible exception of Chern insulators, and 
whenever the time reversal symmetry holds and can be used as a starting point for the inclusion 
of local fields and excitonic effects within time-dependent DFT for such systems. 

In case of metals instead we show that, starting from the density-density correlation function, 
the term which describes the anomalous Hall effect is neglected giving a wrong conductivity. 

PACS numbers: 71.15.-m,71.45.Gm,78.20.Bh 
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I. INTRODUCTION 



The Magneto-Optical Kerr Effect (MOKE) consists in the rotation of the polarization 
plane of light reflected from the surface of a magnetic material. It was discovered in 1877 by 
John Kerr->2 while he was examining the light reflected from a polished electromagnet pole. 
Very recently it became the object of an intense experimental investigations, mainly for two 
reasons. First one can exploit this effect to read suitably magnetically stored information 
using optical means in modern high density data storage technology 3 ' 4 . Second, the MOKE 
can be used as a powerful probe in many fields of research such as microscopy for domain 
observation, surface magnetism, magnetic interlayer coupling in multilayers^ - -. It can also 
be used to observe plasma resonance effects in thin layers, and structural and magnetic 
anisotropics^ - — . 

The microscopic origin of the Kerr effect is a combined action of the spin-orbit coupling 
(SOC) and the net spin-polarization of the material. Indeed the existence of a non zero 
magnetization in the ground state is due to a spontaneous symmetry breaking of the system. 
This symmetry breaking is transferred, through the SOC, to the spatial part of the wave- 
functions so that ?/> + l z (x) is different from ?/>_£ z (x). Accordingly the absorption is different 
for light with right and left circular polarization. 

The problem has been addressed in the literature and ab-initio calculations, based on 
density functional theory (DFT), are available for transition metals like Fe, Co and Nir^-~— 
and, recently, for other materials such as full-Heusler films and Mn-doped GaAs^— . 

The Kerr parameters are commonly obtained from the Kubo formula 3 -^ for the optical 
conductivity tensor using the single particle Kohn-Sham (KS) wavefunctions. This is equiv- 
alent to the computation of the dielectric constant at the random phase approximation 
(RPA) without the inclusion of local fields (LF) and exchange-correlation (xc) corrections. 
We will refer to this approach as independent particles RPA (IP-RPA) scheme. 

Also an alternative approach, based on the Luttinger's Formula^, has been proposed in 

the literature^ 3 - - —, starting from the current-current correlation function, \ ■ Also in this 

=jj 

case the KS wave-functions are used and LF and xc-effects are not considered. 

In all these works the Kerr paramenters are computed within the DFT framework strart- 
ing from the dielectric tensor, e(u), or, which is the same, from the conductivity, g_{uj) 
. However while the diagonal terms of the dielectric tensor, i.e. Eu(u), can be expressed 
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in terms of the density-density correlation function, x° PP ( or Xpp if LF an d xc effects are 
included), at the best of our knowledge, no such expression has been derived for the off- 
diagonal terms, i.e. £ij(uj) with % ^ j. The latter can be obtained only starting from x° 
which however is not expressed as a functional of the density. In this case the sole formal 
justification to the use of KS quantities to construct x° is that "a posteriori" the approach 
give good results and that KS wave-functions can be regarded as a good approximation to 
quasi-particle wave-functions. 

In the present work instead we derive an expression for the full dielectric tensor in terms 
of x°p P an d so for the construction of the Kerr parameters within a density based approach. 
Besides a formal justification to the use of the DFT approach, the result we propose can be 
regarded as a starting point to go beyond the RPA-IP scheme to include LF and xc-effects. 
Indeed while the description of transition metals within the RPA-IP approach is reasonable, 
for semiconductors important deviations are expected, in particular when excitons (magnetic 
excitons) exist. Magnetic semiconductors are in fact materials of great interest, especially 
in view of spin electronics (spintronics) applications and the MOKE can be a valuable tool 
for the investigation of their properties". 

In the present work we limit our discussion to the polar geometry, that is when the prop- 
agation direction of the photon (the z axis) and the magnetization of the system are both 
perpendicular to the surface of the sample (xy). Experimentally this is the most studied ge- 
ometry, and it is also the one which in general gives the largest MOKE signal. To support our 
theoretical derivation with numerical results we have implemented in the plane-waves and 
pseudo-potentials based code Yambo^ the computation of the Kerr parameters. We have 
tested our implementation in transition metals for which well assessed all electrons calcula- 
tions and experimental results are available. For these materials the IP-RPA approximation 
is sufficient. 

Thus in Sec. [Ill we show results on bulk iron, cobalt and nickel in order to validate our 

pseudo-potential based approach, at the IP-RPA level. 

Then we discuss how to construct a formal equation for the dielectric tensor starting 

from x°p P m Sec. IIHI This approach is compared with the one based on x° ■ The two differ 
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by a term which, in general, is zero in systems with an electronic gap or whenever the 
pure time reversal symmetry holds. This term describes the anomalous Hall conductivity. 
Taking iron as reference system we show that, neglecting this term, a large error is induced 
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in the computation of the off diagonal conductivity in metals. However the anomalous 
Hall conductivity is zero in systems with an electronic gap 16 , thus our approach is exact 
for dielectrics. We then show how LF and xc-effects can be included replacing x° pp with 
X pp - The results is a scheme, in principle exact, to compute the Kerr parameters within 
time-dependent DFT (TDDFT). 



II. MOKE PARAMETERS WITH THE IP RPA APPROXIMATION 



A. Theoretical background 



The description of the MOKE can be obtained in terms of the dielectric function e(u) or 
equivalently of the optical conductivity a{u). The two are related by the equation 

to , , 
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where the dielectric tensor at the IP-RPA level can be constructed from the (paramagnetic) 
Xjj, according to the equation^: 
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Here e is the electron charge and m the electron mass, Q the unit cell volume, a labels the 
Cartesian axis, and n = N e i/Q the number of electrons per unit volume. The IP response 
functions is: 
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where the first term on the r.h.s. is the resonant part, while the second in the anti- resonant 
one. c (v) are conduction (valence) band indexes, are the occupation factors, 6j(k) 
the electronic energy at k and Pc„ k (q) the expectation value of the momentum operator 
(ck — q|p a |wk) which in our pseudo-potential based scheme must be computed as^ 



im , 



(ck-q| [p«-—[x\V NL \ 



\vk), 



(4) 



where x a is the a component of the position operator and Vnl the non-local part of the 
pseudo-potential. The infinitesimal rj factor implies that the electromagnetic field is adia- 
batically turned on at t = — oo but can also be viewed as a finite lifetime broadening which 
accounts for scattering process and finite experimental resolution. In the present work we 
use f](u) = 0.3 eV + 0.03hu to mimic an experimental resolution which decrease linearly 



with energy as in Ref. 
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B. MOKE spectra for transition metals 

As a first step we check how a pseudo-potentials based approach preforms in the descrip- 
tion of the Kerr parameters, as it is common wisdom^ that all electrons calculation are 
needed to describe the wave-function in the core region, where the SOC is mostly effective, 
and thus to evaluate the MOKE. 

We start from a ground state DFT calculation for bulk iron, cobalt and nickel with the 
Abinit code^i, using norm-conserving HGH 42 pseudo-potentials and including the SOC; we 
found out that it is crucial, for a correct description of the MOKE, to have the SOC effect 
included both in the pseudo-Hamiltonian and in the construction of the pseudo-potential. 

Bulk iron is studied in its bcc phase with the experimental cell parameter a = 2.87 A; 
an energy cutoff of 65 Ha and a k-points sampling of the Brillouin zone (BZ) 14x14x14. 
Bulk cobalt is studied in its fee phase with the experimental cell parameter a = 3.55 A; 
an energy cutoff of 55 Ha and a k-points sampling of the BZ 8x8x8. Finally bulk nickel is 
studied in its fee phase with the experimental cell parameter a = 3.52 A; an energy cutoff 
of 65 Ha and a k-points sampling of the BZ zone 14x14x14. Semi-core electrons are also 
included in the pseudo-potentials for all systems, as we found the density of states to be 



poorly described using the t 
from the parameters of Ref. 



G 



i pseudo-potentials with only valence electrons, constructed 



42 



Then we compute the dielectric function with the Yambo code^ according to a modified 
version of Eq. [2j 

1 + -rxU(M)J <W - -r xL„(o, w ). (5) 

where the diamagnetic term has been replaced by the zero frequency value of Xjy Indeed in 
cold semi-conductors the diamagnetic term must be exactly balanced by X^j^O, 0), due to 




Figure 1. (color online) Plot of a xx {oj) for bulk BCC iron (Panels a-b), bulk fee cobalt (Panels c-d) 
and bulk fee nickel (Panels a-b). The continuous (red) line are the results from the preset work. 



The dashed line are all electrons resu 
Panels a-b: Ref. 34. Panels c-d: Ref. 
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la . The symbols are experimental measurements. 
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the effective-mass sum rule^ — , but this balance is slowly converging with the number of cv 
states included and the use of xJLj^O) 0) speeds up the convergence^. In metals instead the 
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Figure 2. (color online) Plot of a xy (uj) for bulk BCC iron (Panels a-b), bulk fee cobalt (Panels c-d) 



and bulk fee nickel (Panels a-b) The continuous (red) line are the results from the preset work. The 



dashed line are all electrons results from Ref. 
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difference between the diamagnetic term and the zero frequency value of Xjj gives the Drude 
term, which, with this choice, is set to zero. However it has been shown^ that in practice 
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a ultra-fine sampling of the BZ would be needed to compute this difference. Hence it is 



47. The Drude term 



preferable to include it with a semi-classical model as described in Ref. 
is included only in the computation of the diagonal part of e(u). The optical conductivity is 
finally constructed from Eq. [U Here for the Drude term we used the same parameters of the 
reference all-electrons calculation, i.e. up = (4.9 + zl.87r) eV for iron, up = (8.3 + i27t) eV 
for cobalt and up = (7.5 + i2.24ir) eV for nickel. 

The results for the optical conductivity are plotted in Figs. [T]l2j For all systems there 
is a systematic blue shift of the theoretical peaks against the experimental data. This is a 
known problem of the LDA, due to the self-interaction error, which tends to delocalize the 
(i-orbitals and accordingly gives wrong eigenvalues. The same consideration also explains 
the overestimation of the intensity, as derealization increases the orbitals overlap and thus 
the intensity of the dipoles computed to construct the dielectric function. However a good 
agreement is found with the reference all-electrons calculations. The diagonal component, 
a xx (co), is commonly computed from pseudo-potential based calculations, provided that the 
dipoles are constructed as in Eq. [3l For the off-diagonal component, a xy (u), instead, it 
has been reported that it must be computed using all-electrons wavefunctions^, because 
it depends crucially on the correction to the wave-function due to the SOC term in the 
Hamiltonian. However, in our results, the differences in the diagonal and the off-diagonal 
part of the optical conductivity compared to the reference all-electrons calculations are of 
the same order. Hence we can conclude that pseudo-potentials based calculation can be used 
to compute the Kerr parameters with the same level of confidence of absorption spectra, 
which depend only on the diagonal part of the optical conductivity. We mention that, in the 
literature, pseudo-wavefunctions has been used to construct the off-diagonal conductivity 
at u> = in the computation of the anomalous Hall conductivity^"—. Also in this case a 
good agreement with the all-electrons calculations was found. 

Thus we finally compute the complex Kerr parameters according to the equation 

V K (u) = 9 K {u) + i lK {u) = - — ~ e °» — , (6) 



which is the standard expression for the polar geometry in the small angles limit. Here the 
photon propagates along the z direction and describes a linearly polarized wave with the 
electric field along the x direction. Results are reported in Fig. [3j The blue shift of the 
theoretical results is still present, while the overestimation of the dipoles in both a xx {u) and 
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Figure 3. (color online) Plot of the Kerr paramenters ^k(^) = &K + i^K for bulk BCC iron (Panels 
a-b), bulk fee cobalt (Panels c-d) and bulk fee nickel (Panels a-b). The continuous (red) line are the 
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a xy (u) is compensated and the intensity of the MOKE signal is closer to the experimental 
data than for the case of the optical conductivity. 

To conclude this section, we have shown that for Fe, Co and Ni the Kerr computed 
from our pseudo-potential approach are in good agreement with the results obtained from 
all-electrons calculations. 

III. BEYOND THE IP RPA APPROXIMATION 
A. A density based approach. 

In the previous section we have constructed the Kerr paramenters starting from the KS 
wave-functions using Eq. [2j as it is commonly done in the literature. However the use of 
KS wave-function to construct x° is n °t formally justified. Moreover the inclusion of LF 
and xc effects within a density based approach in Eq. [2] is not straightforward. 

However, to describe the MOKE only the long wave-length term, i.e. q = 0, of the 
dielectric function is needed, where the distinction between longitudinal and transverse 
fields disappears. In this limit the diagonal part of the dielectric tensor can be constructed 
from x° pp which, at finite q, describes only the longitudinal term of the dielectric tensor. 
This approach is formally justified within a density based approach and moreover would 
allow a straightforward inclusion of LF and xc effects within the TDDFT scheme. It is then 
tempting to try to construct the full dielectric tensor at q = from x° pp an d use the result 
to go beyond the IP-RPA scheme. 

Here we provide an heuristic derivation where only longitudinal fields are considered, 
as our final goal is to take the q — » limit. We will prove a posteriori that the result is 
correct for systems with an electronic gap or, more in general, when the pure time-reversal 
symmetry exist and we will discuss in detail the difference between the derived equation and 
Eq. El 

We consider a non uniform system. The dielectric function is defined as: 

W*{q 3 u)=g(qrf,u,)E**(<f,u,). (7) 
Assuming that only longitudinal fields exist, Eq. ([7]) can be written in terms of the potentials 

V ex \q,u) = q|(qq',u;) qY w (q», (8) 
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where the two are related by the equation 
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V**(q, U ) - (q,q>)v to^ (q>)■ 



(9) 
(10) 



Inserting Eq. [10] into Eq. [8] and taking the limit q — > q' — > we can define a generalization 
of the relation that holds between e aa and Xpr^'- 

Aire 2 



lim 



(11) 



In order to compare Eq. [TJJand Eq. [2] we first notice that the latter is divergent for u> — > 0. 
After some algebra Eq. [2] can be rewritten as^ 3 -: 



cvkO 



(12) 



A a/ 3 describes the contribution from the electrons at the Fermi surface, i.e. the Drude 
term, and is zero in cold semiconductors, when there are not partially filled bands. This 
term is also included in Eq. (TT] in the q — >■ limit as discussed in Ref. |47t Once the u 2 



has been isolated using the relation x 
together with 



a 



lip. xJUq.q',") = n p n 5Z / d3k 
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frw - (e ck _ q - e^k) + 277 

(*q ■ xjck)(*q' ■ x ^k) /«k-i(i - /ck) 



(13) 



we obtain the remaining part of Eq. [TTJ 



B. The anomalous Hall effect 

Hence term B a p is not included in Eq. [TTJ It can be explicitly written, at the RPA-IP 
level, as: 



he2 V [d*k(U-f wk ) f^*^ 2 . (14) 



2n 2 m 2 

uw 



This can be shown to be zero when the time-reversal symmetry holds 43 or in any case 
when a = (3 inverting the mute indexes u and w in the second term on the r.h.s. . In 
MOKE experiments however the time-reversal is broken by the existence of a ground state 
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magnetization and by the SOC term in the Hamiltonian and for the construction of ^>k{^) 
we need the terms a ^ (3. 

Thus B a p can differ from zero. In the following, we briefly discuss its physical meaning. 
To fix the ideas we chose a = x and (3 = y. It can be easily proved that the the B xy coefficient 
is (apart a trivial factor raising from the relation between e and a) the intrinsic anomalous 
Hall conductivity (AHC), which is responsible of anomalous Hall effect in magnetic metals. 

In fact, according to Ref.— the AHC reads 

^ - - e jX S Sf unm ' (15) 



i.e. a AHC can be expressed as a BZ integral of the Berry curvature of the u-band, fi^(k) 
(summed over all the occupied states). The latter quantity can be written in terms of the 
ingredients of Eq. [14] as: 56 - 

After some straightforward algebra, one can easily prove that the AHC can be expressed 

as 

Affc B xy 

y =4^' (17) 
which provide the relations between B and the AHC. In the case of magnetic metals, our 
expression constitute an alternative approach to compute a AHC with respect to the methods 
based on the computation of Berry phase- 5 . In the case of insulators instead a AHC has 
been recognized as a topological invariant^, also called Chern number, which can take only 
integer values Thus, in dielectric, B xy can be non-zero only in the so-called Chern insulator, 
hypothetical materials showing a quantum Hall effect without external magnetic field. In 
practice for all the presently known dielectrics eq. [11] can be considered exact. 

Also in this case we have tested, at the IP-RPA level, the effect of B a p on bulk iron com- 
paring the conductivity computed starting either from Eq. [2] or from Eq. [TT] In Fig. [H (a — b) 
we show the error induced in the computation of the off-diagonal conductance on bulk iron. 

To clarify the relation between this difference and the AHC also numerically we have 
considered, in Fig.[U(c — d), the plot of the conductivity at small smearing, i.e. r] = 0.05 eV, 
as Eqj2] and Eq. [TTJ are equal only in the limit 77 — ^ — . From Fig. [H (c) is clear that 
the difference of the two gives a constant value, as expected theoretically, a part from 
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Figure 4. (color online) Off-diagonal element of the conductivity tensor, cr xy (u), of iron computed 
starting from Eq. [Til dot-dashed (green) line and from Eq. [2] continuous (red) line. The dashed line 



are all electrons results from Ref. 
Ref. 



35l ; empty squares, Ref. 



la . The symbols are experimental measurements: filled circles, 
36l . In panels a-b a xy {oj) is computed with the smearing used in the 
present work. In panels c and d instead a xy {oj) is evaluated at rj = 0.05 eV. The difference between 
the dot-dashed (green) line and the continuous (red) line is represented with a thin (black) line. 



In the inset of the panels also the difference in terms of the dielectric function is represented. 



the region oj ~ 0, where the 1/u 2 term makes Eq. [2] numerically unstable. We can thus 
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extract the cr^ c = 665 (Qcm) 1 which is not so far from the theoretically computed value 



750.8 (Qcm)~ l of Ref. [56|. The difference is likely due to the sampling of the BZ. As for 
the case of the Drude term, the anomalous Hall conductivity depends on the contribution 
from the electrons at the Fermi surface and thus a very fine sampling of the BZ should be 
needed, which is beyond the scope of the present work. Also we see in Fig. Hl(cf) that at 
small 77 the difference between the imaginary parts of the conductivity computed with the 
two approaches goes to zero (Fig. HI (6)) as expected from the theoretical derivation. Finally 
in the insets we have also represented the difference at the level of the dielectric function 
which are Im[Ae] ocl/u and Re[Ae] oc 8(u) thus respecting the Krames-Kronig relations. 

C. Inclusion of local fields and excitonic effects. 

The generalization of Eq. [2] to include LF and xc-effects is non trivial and needs a 

careful distinction between longitudinal and transverse induced fields. The result, derived 

in Ref. 59(, is to replace the IP %° with the one constructed from the analytical part of the 

_=jj 

electron-hole (eh) propagator L(12, 34) solution of the modified Bethe-Salpeter equation: 

1(12, 34) = L (12, 34) + L (12, 1'2') 

v(l', 305(1', 2')8(3', 4') - iW(l, 2)5(1', 3')8(2', 4')j L(3'4', 34), (18) 

with 1 representing spatial, time and spin coordinates: 1 = (xi,ti,o"i). The long range part 
of the exchange interaction v(l, 2) = 5(t 1 — t 2 ) / l x i — x 2| between the electron and the hole 
is truncated with the substitution v — > v where in reciprocal space 

47re 2 /|q+G| 2 if G ^ 0, 

Mq)H ■ (is) 

if G = 0. 

From the electron-hole propagator the \ and x are constructed with the relations^ 2 - 

X pp (l,2) = -z^L(l,2;l + ,2 + ), (20) 
|..(1,2) = -ih^ [(Vi - Vi)(V a - V' 2 )I(1,2; l'2')] 1/=1+ 2 , =2+ , (21) 

with 1 + = /im T _ i .o( x i, ti + r, a±) The result is then 

ScA") = [ 1 - -^jr ) 5 ^ - — X Jaj0 (O,u). (22) 
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A possible strategy to remain within a density based formalism—, starting from Eq. [2] 
could be to use the unphysical L TDDFT replacing iW(l, 2)5(1', 3')6(2', 4') with f xc (l, 2)5(1', 2')6(3 
in Eq. [TSJ However this is not formally justified neither and at least a current based for- 
malism should be used, i.e. current-DFT— ~— . Indeed for the description of absorption 
spectra, the diagonal part only of the dielectric function is commonly constructed from \ pp 
to include LF and xc effects^ starting from £u(uj)[x pp }- A similar derivation can be used 
to include LF and xc-effects replacing x° pp with \ pp in Eq. [TTJ In this case however the 
Dyson equation for the response function should be written for a non homogeneous system 
assuming, as we did in the IP case, that transverse fields can be neglected as we are looking 
for the q — > limit. The result is: 

Aire 2 

£ a p(u) = 5 aP - lim ——x(q a , qp, w), (23) 
q a ,g^o q z " 

which, according to the discussion of the previous sections should hold when the time- 
reversal symmetry exists or for systems with an electronic gap 16,65 . Eq. [23] must then be 
compared with Eq. [22] If x and x are constructed from the same L using Eqs. I20H2T1 
the two are diagonalized by the same vectors in c v space, A 1 ^; here I is the index of the 
excitation, which now can be a mixture of electron-hole pairs. In this case inserting the 
vectors in the equations the two approach will differ by the term 

BaP = 2^ 1* 2^ J d k(/ " k " U) {h^f ' (24) 

I uw J 

which defines a generalization of the Anomalous Hall effect. Here Uj are the poles of X pp - 
In common metals usually we have A ! cvk = 5/,(cu)i (i-e. each vector A 1 is different from zero 
only for a specific transition cv = (cv)i ) and ftwi = e c k — e v k, thus Eq. EH reduces to Eq. [T4l 
However if one remains within a pure DFT approach, then the vectors which diagonalize 
Xppi ^Ivk 00 ^ do not, in general, diagonalize L and thus % . In this case Eq. [23] and Eq. [22] 
could also differ by a term proportional to A^ — A l cv ^. This term must be zero for 
a = (3, while its relevance in the case a ^ (3 and its eventual physical meaning are left under 
study. 



IV. CONCLUSIONS 



We have proposed a scheme to compute the magneto-optical Kerr effect in magnetic- 
semiconductors. The scheme has two main novelties. First is based on pseudo-potentials 
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calculations. This is the most widely used approach to describe extended systems and we 
have shown that pseudo-wavefunctions can be used to obtain the Kerr parameters. The 
results we find are comparable with all-electrons calculations, provided that the Spin-Orbit 
interaction is correctly accounted for in the construction of the pseudo-potential. 

Second we have discussed the inclusion of local-field and excitonic effects in the computa- 
tion of the MOKE. We have shown that two strategies can be used: (i) the Bethe-Salpeter 



equation, through the result derived in Ref. |59|, but also, in almost any case of interest, 



(ii) an approach based on time-dependent density-functional theory and in general on the 
density-density correlation function through the result derived in the present manuscript. 
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